Effective Polyakov Loop Dynamics for Finite Temperature G 2 Gluodynamics 
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Based on the strong coupling expansion we obtain effective 3-dimensional models for the Polyakov loop in 
finite-temperature Gi gluodynamics. The Svetitsky-Jaffe conjecture relates the resulting continuous spin mod- 
els with G2 gluodynamics near phase transition points. In the present work we analyse the effective theory 
in leading order with the help of a generalised mean field approximation and with detailed Monte-Carlo sim- 
ulations. In addition we derive a Potts-type discrete spin model by restricting the characters of the Polyakov 
loops to the three extremal points of the fundamental domain of Gi. Both the continuous and discrete effective 
models show a rich phase structure with a ferromagnetic, symmetric and several anti-ferromagnetic phases. The 
phase diagram contains first and second order transition lines and tricritical points. The modified mean field 
predictions compare very well with the results of our simulations. 
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I. INTRODUCTION 

For many good reasons G-2 gluodynamics has attracted much 
attention recently. For example, the 14-dimensional excep- 
tional group G2 has a trivial centre, in contrast to the usually 
studied SU(N) gauge groups. Indeed it is the smallest simple 
and simply connected compact Lie group with this property. 
Thus Gi gluodynamics is useful to better understand the rel- 
evance of the centre symmetry for confinement [ 1 ] . Actually 
a non-trivial centre is needed in several proposed scenarios 
for confinement and hence G2 gluodynamics can be used to 
test these proposals. It has been convincingly demonstrated 
that the theory shows a first order finite temperature transi- 
tion without order parameter from a confining to a deconfin- 
ing phase which can be explained by centre vortices [2]. In 
this context confinement refers to confinement at intermedi- 
ate scales, where a Casimir scaling of string tensions has been 
reported [3]. But on large scales, deep in the infrared, strings 
break due to gluon production and the static inter-quark poten- 
tial becomes flat [4]. Recently it has been demonstrated that 
chiral symmetry is broken at low temperatures and is restored 
at high temperatures at the thermodynamic phase transition 
[5]. 

G2 gluodynamics has an intriguing connection to SU (3) 
gauge theory. When one couples a scalar field in the 7- 
dimensional fundamental representation to the gauge field 
one can break the G2 gauge symmetry to the SC/(3) gauge 
symmetry of strong interaction. With increasing hopping pa- 
rameter k the resulting Yang-Mills-Higgs theory interpolates 
smoothly between G2 gluodynamics without centre symmetry 
and SU(3) gluodynamics with Z3 centre symmetry. For in- 
termediate values of the hopping parameter the theory mimics 
ST/ (3) gauge theory with dynamical quarks and the masses 
of these 'quarks' increase with increasing hopping parameter. 
In G2 gluodynamics the Polyakov loop is no longer an order 
parameter in the strict sense. Despite this fact it still serves as 
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FIG. 1 : Expectation value of the traced Polyakov loop in the funda- 
mental 7-dimensional representation in G2 gluodynamics on a 16 3 x6 
lattice as obtained via hybrid Monte-Carlo sampling. The phase tran- 
sition is located for the Wilson action at a critical j3 c ~ 9.765. 



an approximate order parameter separating the confined from 
the deconfined phase (see Fig. 1) with a rapid change at the 
phase transition point. 

According to the conjecture by Svetitsky and Yaffe [6, 7] 
the dynamics at the finite temperature confinement-decon- 
finement transition of a d+ 1 -dimensional pure gauge theory 
can be described by an effective spin model in d dimensions. 
Based on our earlier results on finite-temperature SU (2) and 
ST/(3) gluodynamics [8-12] there are strong indications that 
the effective models derived and analysed in the present work 
are sufficient to accurately describe the dynamics of Polyakov 
loops. The direct connection between the effective spin mod- 
els and G2 gluodynamics is postponed to a forthcoming pub- 
lication. 

In Sec. II we review kinematic aspects of G2 and the main 
implications for G2 gluodynamics. Afterwards in Sec. Ill the 
strong coupling expansion for the effective Polyakov loop ac- 
tion is explained and in particular the effective theory in lead- 
ing order is introduced. In Sec. IV we investigate the proper- 
ties of the effective model first by a classical analysis, then by 
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a modified mean field approximation and finally by extensive 
Monte-Carlo simulations. Reducing the continuous spin de- 
grees further to the discrete spins situated at the 3 edges of the 
fundamental domain of G2 we end up with a deformed Potts- 
type spin model whose phase diagram is explored in Sec. V. 



H. THE GROUP G 2 

G-2 is the smallest of the five exceptional simple Lie groups 
and can be viewed as a subgroup of SO (7) subject to seven 
independent cubic constraints for the 7-dimensional matrices 
g representing SO (7) [13]: 

T a bc = Tdef 9da 9eb 9fc- (1) 

Here T is a total antisymmetric tensor given by 

T\27 = T154 = T163 = T235 = ^264 = T374 = Ts7g = 1. 

(2) 

The constraints (1) for the group elements reduce the 21 
generators of SO(7) to 14 generators of the group G2 with 
rank 2. Its fundamental representations are the defining 7- 
dimensional and the adjoint 14-dimensional representation 
with Dynkin labels 

(7) = [1,0], (14) = [0,1]. (3) 

G-2. has a trivial centre and its Weyl group is the dihedral 
group Dq of order 12. Additionally G2 is connected to SU (3) 
through the embedding of SU(3) as a subgroup of G2 accord- 
ing to [14] 

G 2 /SU(3) ~ SO(7)/SO(6) ~ S 6 . (4) 

So when the 56 part of G2 is frozen out 1 we end up at SU (3) 
gauge theory. 

In effective theories for the gauge invariant (traced) Polya- 
kov loops in the fundamental representations we are aiming 
at, only the reduced Haar measure is needed. Based on [15] 
this measure can be given for a parametrisation of the conju- 
gacy classes either by angular variables or alternatively by the 
fundamental characters, 

djjL oc J 2 d(p\ d(f2 — J d\7 dxiA- (5) 

The density J 2 can be expressed in terms of the fundamental 
characters, 

J 2 = (4x? - X? - 2X7 - 10X7X14 + 7 - 10xi4 - xL) 
x (7 - x? - 2X7 + 4xi 4 ) , 

(6) 



This is possible when a fundamental Higgs field is coupled to the gauge 
field [1]. 



where the characters are given in terms of (particularly cho- 
sen) angular variables ip\ t2 as 

X 7 = 1 + 2cos(v?!) + 2cos(^ 2 ) + 2cos(<£i + ip 2 ), 

X14 = 2(l + COs(l^l) + C0S(<^1 - tp 2 ) + COs(lf 2 ) 

+ cos(<^i + ip 2 ) + cos(2<ySi + ip 2 ) + cos(<pi + 2tp 2 )). 

(7) 

The boundary of the fundamental domain is determined by 
J = and thus is parametrised by the three curves (see Fig. 2) 



Xi4 = -5(x7 + l)±2(x7 + 2) 3 / 2 . 

Note that the reduced G2-Haar measure is maximal not at the 
origin but for (x7, X14) = (~ 1/5,-2/5). The fundamental 
domain has no symmetries at all and this expresses the fact 
that the centre of G 2 is trivial. 




FIG. 2: Fundamental domain of G2. Darker regions indicate a bigger 
Haar measure. 



Representation theory and implications for confinement 

In the pioneering work [13] the confining properties of G2 
have been discussed and compared to those of SU (3). Quarks 
and anti-quarks in SU(3) transform under the fundamental 
representations 3 and 3 such that their charges can only be 
screened by particles with non-vanishing 3-ality, especially 
not by gluons. This explains why in the confining phase of 
SU(3) gluodynamics the static inter-quark potential is lin- 
early rising up to arbitrary long distances. As a consequence 
the free energy of a single quark gets infinite and the Polyakov 
loop expectation value vanishes. Hence the Polyakov loop 
discriminates the confining from the deconfining phase and 
at the same time serves as order parameter for the Z3 centre 
symmetry. 

To better understand G2 gluodynamics we recall the de- 
composition of tensor products into irreducible representa- 
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tions, 

(7) <g> (7) = (1) © (7) © (14) © (27) 
(7) $ (7) © (7) = (1) © 4 • (7) © 2 ■ (14) © 3 ■ (27) 
© 2 • (64) © (77) 
(14) © (14) = (1) © (14) © (27) © (77) © (77') 
(14) © (14) © (14) = (1) © (7) © 5 • (14) © 3 • (27) i 



with Dynkin labels 



(1) 
(77) 



[0,0], 
[3,0], 



(27) 
(77') 



[2,0], 
[0,2]. 



(64) = [1,1], 



(9) 



(10) 



The quarks in G2 transform under the 7-dimensional funda- 
mental representation, gluons under the 14-dimensional fun- 
damental (and at the same time adjoint) representation. From 
(9) we see that similarly as in SU(3) two or three quarks can 
build a colour singlet (meson or baryon, respectively). In G2 
gluodynamics three centre-blind dynamical gluons can screen 
the colour charge of a single quark, 



(7) © (14) © (14) © (14) = (1) 



(ID 



Thus the flux tube between two static quarks can break due to 
gluon production and the Polyakov loop does not vanish even 
in the confining phase [1]. This shows that the Polyakov loop 
can at best be an approximate order parameter (see Fig. 1) 
which changes rapidly at the phase transition and is small (but 
non-zero) in the confining phase. To characterise confinement 
we can no longer refer to a non-vanishing asymptotic string 
tension and vanishing Polyakov loop. Instead we define con- 
finement as the absence of free colour charges in the physi- 
cal spectrum. In the confining phase the inter-quark potential 
rises linearly at intermediate scales [2, 3]. 



III. EFFECTIVE THEORIES AND THE STRONG 
COUPLING EXPANSION 

Based on a conjecture relating finite temperature SU (N) glu- 
odynamics in d+1 dimensions at the critical point with a Zn 
spin model in d dimensions [6, 7], there have been extended 
studies to compare correlation functions of both systems for 
SU{2) [9-11] and SU{3) gluodynamics [8, 12], either by us- 
ing Schwinger-Dyson equations or demon methods [16, 17]. 
The strong coupling expansion for the distribution of the in- 
homogeneous Polyakov loops was taken as ansatz for the (ex- 
ponentiated) effective Polyakov-loop action. This way effec- 
tive models for SU (3) gluodynamics have been derived in [8]. 
Here we sketch how one arrives at the analogous results for 
G-2 and obtain the effective continuous spin model in leading 
order. 

Starting with the lattice Wilson action 



= 



6 



(12) 



a strong coupling expansion (for small (J) is performed to ar- 
rive at an effective theory for the local Polyakov loops. To 



do that one inserts a group valued delta function into the path 
integral, 

Z = J T>Uexp(-Sw[U}) 

= J VP J VU 5 (p»,JlU T ^ cxp(-S w [[/]) (13) 

= J W exp(-5 e ff['P]). 

Here VP denotes the product of reduced Haar measures on 
the sites of the spatial lattice. We do not need the full Haar 
measure of G2 since the effective action S e s only depends on 
the gauge invariant content of the local Polyakov loop. 

In compact form the strong coupling expansion is then 
given by 



r Ki...K T ti...ir i=l 

i 

with the basic building blocks 



(14) 



\x-y\=l. (15) 



Here r counts the number of link operators contributing at 
each order. The coefficients c^ lT n couple the operators 
Snuii sitting at sites Xi and yi separated by distance ti in 
representation TZi. The effective action hence describes a net- 
work of link operators that are collected into (possibly discon- 
nected) 'polymers' contributing with 'weight' c^" ^ ■ The 
resulting 'operators' (Polyakov loop monomials) are dimen- 
sionless and there is no natural ordering scheme at hand. Our 
chosen truncation is based on the strong coupling expansion in 
powers of j3 which is closely related to the dimension of the 
corresponding group representations and the distance across 
which the Polyakov loops are coupled. In the strong coupling 
expansion truncated at 0{(3 kNt ) one has r < k and the addi- 
tional restriction \TZi | +■ ■ ■+ 1 lZ r | < k with \R\ =pi +P2 for a 
given representation 1Z of G2 with Dynkin labels [pi , P2]. The 
leading order terms only contain interactions between nearest 
neighbours (xy) and the two fundamental representations. 

For SU(3) the characters of the two fundamental represen- 
tations are complex conjugate of each other such that the ef- 
fective Polyakov loop action contains just one term in lead- 
ing order. In G2 this situation changes and we find two inde- 
pendent contributions in leading order. We refer to the corre- 
sponding model containing the two fundamental representa- 
tions as fundamental model. Its action is explicitly given by 

S eS = MJ2 Xi{V x )x7{V v ) +A 14 J2 XM(V x )xi4(Vy), 

(xy) (xy) 

1 ■> s ' 

SV S14 

(16) 

where the couplings A 7 and A14 are indexed by the dimension 
of the involved representation. In next-to leading order there 
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exist six additional terms with nearest neighbour interactions. 
Their explicit forms are dictated by representation theory (9) 



S27 = Xn{V x )X27(V y ), 
(ay) 

5*77' = Y X77'CP x )x77'CPy), 
#64 = Y XGii'PxjXSii'Py), 

{ny) 



(17) 



(xy) 



514,14 = Y (Xl4(V x )xi4(Vy)) 2 , 
(a; y) 

S 7 ,U = Y, X7{V x )X7{Vy)Xli{V X )Xli(Py). 

(ny) 

It the remainder of this work we shall neglect the next-lo lead- 
ing order terms and concentrate on the fundamental model 
(16). 

IV. THE FUNDAMENTAL MODEL 

For the fundamental effective model (16) we shall localise 
the symmetric, ferromagnetic and anti-ferromagnetic phases 
with coexistence lines in order to find the region in the space 
of couplings A7, A14 where a connection to G2 gluodynamics 
can be established. 



A. Classical analysis 

For strong couplings the fluctuations of the Polyakov loops 
are suppressed and the spin system behaves almost classically. 
Thus for large A7I and | A14 1 we may compute the phase di- 
agram by minimising the classical action. Anticipating that 
there are anti-ferromagnetic phases we introduce the odd and 
even sublattices 



A G = {x I xi + X2 + X3 odd} and 
A e = {x I xi + x 2 + x 3 even} . 



(18) 



On each sublattice the Polyakov loop is assumed to have a 
constant value and the two values are denoted by V a and V t , 
respectively. We denote the corresponding characters in the 
fundamental domain of G2 by 



( 



Xo 



X7,e > 


_ ( X7 


xi4, e y 


\Xu 




_ f X7 


xi4, y 


~ Ul4 



0P e ) and 

(■Pa)- 



(19) 



With this assumption and notation the action of the fundamen- 
tal model (16) reads 

SWf = V X jK Xo , with K = 3 A M . (20) 



To localise the different phases we may assume that the 
Polyakov loop on one sublattice, say A Q , is equal to the group- 
identity with maximal characters, X7,o = 7 and xi4,o = 14. 
For given couplings A7 and A14 the corresponding thermody- 
namic phase is then determined by that Polyakov loop on A e 
for which the linear function 7A7X7.e + 14Ai4Xi4 je is minimal. 
If the minimising characters are the same on both sublattices 
then the phase is ferromagnetic, else it is anti-ferromagnetic. 
The minimum of the linear function is attained for Xe on one 
of the three corners of the fundamental domain in Fig. 2 or 
on the curve connecting the corners (—1, — 2) T and (7, 14) T . 
Depending on the sign of A7 and the slope £ = A14/A7 one 
finds the following phases: 

• For A 7 > and £ < -1/2 or for A 7 < and £ > -1/8 
we find the ferromagnetic phase F with xj = (7, 14). 

• For A7 > and £ > 1/14 we find a anti-ferromagnetic 
phase AF1 with xj = (~ 1; _ 2). 

• For A7 > and —1/2 < £ < 1/14 we find a second 
anti-ferromagnetic phase AF2 with xj = ( — 2, 5). 

• For A7 < and £ < —1/8 the characters xj = (~ 1 ~ 
l/£, l/(2£) 2 - 2) change continuously from (7, 14) to 
( — 1,-2) along the connecting boundary curve of the 
fundamental domain. This transition phase is denoted 
by F — > AF1. 

The phase portrait is depicted in Fig. 3 where we also included 
the expected symmetric phase for weak couplings. Since in a 
symmetric phase entropy wins over energy it cannot be seen 
in any classical analysis. 




FIG. 3: The classical phase diagram of the fundamental effective 
G2 model. In addition to the calculated ferromagnetic and anti- 
ferromagnetic phases we expect a symmetric phase for weak cou- 
plings. 



B. Mean field analysis 

The classical analysis is refined by a modified mean field ap- 
proximation in which one allows for inhomogeneous mean 
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fields. First we recall the main aspects of the method. Here 
we are interested in expectation values of the form 



(A) 



vr = \{d^v x ), 



(21) 



where the gauge invariant function A(V) depends on the 
Polyakov loop via the the fundamental characters X7^P) an d 
Xui'P) and d^i is the reduced Haar measure of G2. The equi- 
librium probability measure W cxp(— S e ff) / Z is the unique 
minimum to the variational problem 



inf (S e ff 
p 



lnp) 



(22) 



where the p-indexed expectation value is calculated with the 
integration measure p[V] VP, whose probability density p is 
to be varied. Expectation values of observables can then be 
computed as 



(A) p 



VPp[P]A[P]. 



(23) 



In this scope a Monte-Carlo simulation is just the approxima- 
tion of the probability density p[P] oc exp(—S e g[V]) with 
a finite set of configurations which give p[P] w N, 



MC 



J2t=i ^{"P ~ %)> where Vt is the configuration in the t'th 
MC step and Amc is the number of MC steps. 

In a variational approach the mean field approximation 
amounts to the restriction of the admissible densities p to 
product form 



(24) 



Then expectation values factorise and the computation can be 
done site by site. Due to the translational invariance of the 
action one may believe, that the minimising density is trans- 
lational invariant, p x (T) = p(V). However, this assumption 
is only justified for the symmetric and ferromagnetic phases 
with constant mean fields. 

Anticipating the existence of additional anti-ferromagnetic 
phases we partition the lattice into its even and odd sub- 
lattices, as we did in the classical analysis, and allow for dif- 
ferent densities on the sublattices, 



p x (P x ) = 



PeCPx) : x e A e 
Po{V x ) : x G A G . 



(25) 



The classical analysis is then recovered by allowing only 5- 
type point-measures forp e ,o- In the modified mean-field anal- 
ysis we allow for all p e o in the variational principle with pre- 
scribed mean fields Xe an d Xo on the even an d odd sublattices. 
The effective potential u is then obtained by computing 



"(Xe,Xo) 



— inf (Ses 
V p 



(26) 



subject to the following four constraints for the admitted den- 
sities p e and p a : 



(X>e 



(27) 



The one-site expectation values (. . .) are calculated with 
p e ,o d/i. To actually compute the minimising densities one 
needs the expectation value of the action and entropy, given 
by 



(S ef f) p = VxjKxo 
V 



(In 35) 



V 



( ln Po) : 



(28) 



where K is the matrix given in (20). On each sublattice the 
variational problem is solved by a density p oc exp(j • xi'P)) 
with two Lagrangian multipliers j = (jV, Ju). The four mul- 
tipliers are determined by the four constraints in (27). Using 
this solution for the densities in the variational principle de- 
termines the effective potential as function of the prescribed 
mean fields. The expectation values of the characters on the 
two sublattices minimise the effective potential. These mini- 
mas solve the following system of coupled gap equations 



dw{x ) 



KXo 



dw(Xe 



dXo dx e 
w(x) = ln / dn(V)e-* TKx{v) . 



(29) 



We have calculated the expectation values of X7 an d Xi4 on 
both sublattices as functions of the couplings on a 120 x 100 
grid in the rectangle 

-0.3 < A 7 < 0.3 and - 0.25 < A14 < 0.25. (30) 



The contour plot of the expectation value 



1 



(xt) = ((X7,e) + (X7,o» ■ 



(31) 



called magnetisation, is depicted in Fig. 4. As expected, for 
weak couplings we find a symmetric phase with vanishing 
magnetisation in the centre of the phase diagram. On the 
lower left, for negative couplings, we find the ferromagnetic 
phase with (X7,e) = (X7,o) ~ 7 or equivalently with a typical 
V x near the identity. For an unambiguous identification of the 
phases one needs the expectation values of both xt an d Xi4 on 
both sublattices. We have calculated these four expectation 
values for the fundamental model on a grid in the space of 
coupling constants with extensive Monte Carlo simulations. 
Since the numerical simulations and mean field approxima- 
tion yield almost identical results we defer the detailed dis- 
cussion of phase portrait, and in particular the localisation of 
the various anti-ferromagnetic phases for positive couplings, 
to the following section. 



C. Monte-Carlo results 

We performed our Monte-Carlo simulations with about 10 000 
samples for every point on a 60 x 50 grid inside the rectangle 
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FIG. 4: Polyakov loop (^7) of the fundamental G2 spin model in 
mean field approximation. 



FIG. 6: Staggered magnetisation S of the fundamental model ob- 
tained via MC simulation on an 8 li lattice. 



(30) in the space of coupling constants. Two neighbouring 
points on this grid are separated by 0.01. First we calculated 
the magnetisation (xr) and the resulting phase portrait is de- 
picted in Fig. 5. It looks very similar to the portrait calcu- 
lated in the mean field approximation, see Fig. 4. For weak 




FIG. 5: Polyakov loop (x?) of the fundamental model obtained via 
MC simulation on an 8' 5 lattice. 

couplings entropy wins over energy and the product of Haar 
measures of the Polyakov loops become relevant. In order 
to unambiguously identify the anti-ferromagnetic phases for 
positive couplings we again subdivided the lattice in the even 
and odd sublattice, A = A e U A D , and measured the staggered 
magnetisation 

S = ^<IX7,e-X7,o|>. (32) 

The corresponding contour plot is exhibited in Fig. 6. On 
top and on the right of the plot the staggered magnetisation 
gets large and we identify this region as belonging to anti- 
ferromagnetic phases. For large absolute values of A7, A14 
action (energy) dominates over entropy and this explains why 



the simulation results agree well with the classical analy- 
sis for strong couplings: all phases but the transition phase 
F — > AFl are already visible in the classical phase diagram in 
Fig. 3. In the Monte-Carlo simulation an additional "symmet- 
ric phase" with vanishing Polyakov loop and vanishing stag- 
gered magnetisation appears for weak couplings, in complete 
agreement with our mean field analysis. The resulting phase 
diagram with one symmetric, one ferromagnetic and two anti- 
ferromagnetic phases is depicted in Fig. 7. 
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-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



A7 

FIG. 7: Phase diagram of the fundamental effective model. The or- 
ders of phase transitions are indicated with full lines (first order) and 
dotted lines (second order/crossover). 

Eventually the finite temperature phase transition in G2 
gluodynamics will correspond to a transition between the 
symmetric and the ferromagnetic phase in the effective spin 
model. The dependence of the effective couplings A7, A14 on 
the Yang-Mills coupling g can be calculated with the help 
of powerful inverse Monte-Carlo techniques [8, 12]. This 
will be done in a forthcoming publication. However, we 
anticipate that the confinement-deconfinement phase transi- 
tion in G2 gluodynamics will happen near the critical point 



0.0001 I 1 1 1 1 1 1 1 1 

-0.16 -0.14 -0.12 -0.10 -0.08 -0.06 -0.04 -0.02 0.00 

At 

FIG. 8: Polyakov loop (xi) of the fundamental effective model with 
coupling A14 = obtained via MC simulation on an 8 3 lattice. 



Ai4, c = 0, A7. c ~ —0.0975(75) of the fundamental model. 
Thus we have plotted the magnetisation in the vicinity of this 
first order transition from the ferromagnetic to the symmetric 
phase in Fig. 8. Note that even in the "symmetric" phase we 
find a non-zero magnetisation (^7) which jumps at the critical 
coupling X 7 c . This parallels the jump of the Polyakov loop in 
G2 gluodynamics, see Fig. 1. 

The phase diagram in Fig. 7 contains lines of second and 
first order transitions and 3 triple points. The full lines be- 
long to first order and the dotted lines to second order tran- 
sitions. Note that we may pass from the symmetric to the 
ferromagnetic phase via a first or via a second order tran- 
sition. The transitions from the ferromagnetic to the anti- 
ferromagnetic phases AF1 and AF2 and between the anti- 
ferromagnetic phases are always of second order. In order 
to determine the orders of the transitions we calculated more 
than 30 histograms for the Polyakov loop distribution near 
the various phase transition curves and the changes of vari- 
ous 'order parameters' when one crosses the transition lines. 
A typical scatter plot is depicted in Fig. 9. It shows the dis- 
tribution of X7 at a transition from the symmetric to the anti- 
ferromagnetic phase AF1 with critical couplings A7. c = and 
Ai4, c = 0.1446. Without further analysis it is already clear 
that we are dealing with a first order phase transition. 

The following Fig. 10 shows the behaviour of the magneti- 
sations (xt) an d (X14) near the transition from the symmetric 
to the ferromagnetic phase, which happens for A14 = 0.13 
and A7 between —0.18 and —0.12. Both expectation values 
vary continuously during the transition and this already sug- 
gests that the transition is of second order. This conclusion is 
further substantiated by the corresponding histograms for the 
distributions of \7 an d Xi4- 



V. THE G 2 POTTS MODEL 

After having collected sufficient information to reconstruct 
the full phase diagram of the fundamental continuous spin 



-0.5 1 1 1 1 1 1 1 
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X7 

FIG. 9: Distribution of \7 and Xi4 m me fundamental domain of G2 
at A 7 = 0andA 14 = 0.1446. 
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FIG. 10: Magnetisations {X7) (upper curve) and (xn) (lower curve) 
for various A7 at A14 = 0.13. 



model with two effective couplings we now truncate the de- 
grees of freedom further to arrive at a discrete spin model. In 
the case of the well-studied SU(3) Polyakov loop models one 
projects a Polyakov loop V x onto the closest centre elements 
of SU(3) and arrives at a Z3 Potts model with action (energy) 
given by 

Sn = -0 J~] S(a sc , ay), <7 X G Ztv- (33) 

The continuous and discrete models have coinciding critical 
exponents at the second order anti-ferromagnetic phase tran- 
sition and similar phase structures [18]. Motivated by these 
earlier successes we perform a similar reduction of the funda- 
mental G2 spin model and arrive at a discrete Potts-like G2 
spin model. 

By projecting the values of V x to the three group elements 
with characters X7, X14 lying at the extremal points of the fun- 
damental domain in Fig. 2, we arrive at a model for the tree 
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FIG. 11: The classical phase diagram of the discrete G2 Potts model, 
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with nearest neighbour interaction determined by the Potts- 
type action 

Spotts = a ^ [0 AiJ Gy ■ (35) 

(<»,») 

As expected, the classical phase diagram of the Potts-type 
model with discrete spins is similar to the diagram of the fun- 
damental model with continuous spins. Depending on the sign 
of A7 and the slope £ = A14/A7 we find the following phases 
and phase transition lines: 

• For A7 > and £ > -1/2 and for A7 < and £ > 
—47/206 we find the ferromagnetic phase F with Xe = 
Xo = (7, 14) T . 

• For A 7 < and -1/2 < £ < -47/206 we find the 
anti-ferromagnetic phase AF3 with Xe — ( — 1> — 2) T 
and X o = (-2,5) T 

• For A7 > and -1/2 < £ < 1/14 we find the anti- 
ferromagnetic phase AF2 wit Xe = (~ 2, 5) T and Xo = 
(7, 14) T . 

• For A7 < and £ < -1/2 and for A7 > and £ > 
1/14 we find the anti-ferromagnetic phase AF1 with 
Xe = (-l,-2) T and Xo = (7,14) T . 

The phase portrait is depicted in Fig. 11, where we have 
inserted by hand the expected symmetric phase for weak 
couplings. A striking difference between the diagrams in 
Fig. 3 and in Fig. 1 1 is the absence of the "transition phase" 
F — ► AF1 in the discrete model for which this phase does 
not exist by construction. Instead we find a third anti- 
ferromagnetic phase denoted by AF3 in Fig. 11. In addi- 
tion, in the symmetric phase of the continuous spin model 




FIG. 12: Magnetisation (xr) of the discrete G2 Potts model in mean 
field approximation. 



(X7) ~ and in the symmetric phase of the discrete spin 
model (x 7 ) ~ 4/3. 

Similarly as for the continuous model we calculated the 
phase diagram of the discrete model with the help of the mod- 
ified mean field approximation. The contour plot for the mag- 
netisation is depicted in Fig. 12. In the lower part of the plot 
we can see the ferromagnetic phase for which the Polyakov 
loops on both sublattices A e and A Q are equal to the identity 
with very high probability. 

The corresponding contour plot as obtained from Monte- 
Carlo simulations is shown in Fig. 13. Again, the mean field 
approximation and the Monte-Carlo simulations fully agree 
over the whole range of coupling constants. Note that the clas- 
sical behaviour as depicted in Fig. 1 1 can be seen already for 
rather small coupling constants. 




FIG. 13: Magnetisation (xi) of the discrete G2 Potts model obtained 
via MC simulation on an 8 3 lattice. 

In order to localize the anti-ferromagnetic phases we also 
measured the staggered magnetisation introduced in (32). The 
resulting values on a grid in coupling constant space are plot- 
ted in Fig. 14. In accordance with the classical analysis we de- 
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FIG. 14: Staggered magnetisation S of the discrete G2 Potts model 
obtained via MC simulation on an 8 3 lattice. 



tect different anti-ferromagnetic phases for positive coupling 
constants in the upper right part of the phase portrait. 

As concerning the relation between the Potts-type model 
and G2 gluodynamics one caveat should be mentioned. In the 
discrete spin model there exists no real "symmetric phase" 
with a fixed expectation value of X- Even for very weak 
coupling do the magnetisations (x) depend on the couplings. 
This is a remnant of the missing centre symmetry of G2. Nev- 
ertheless, there exists a first order phase transition from one 
(would be symmetric) ferromagnetic phase at A7 w A14 « 
to a second ferromagnetic phase with V x directed to the 1- 
element in G2. Even in the simple discrete model we see 
very pronounced what happens in real G2 gluodynamics. In 
the deconfining and confining phase there is a non-vanishing 
Polyakov loop, which still serves as approximate order pa- 
rameter for confinement since it shows a steep jump at the 
transition point. 



coupling expansion of G2 gluodynamics. Already the leading 
order continuous and discrete effective theories show a rich 
phase structure with two coexisting phases along transition 
lines and three coexisting phases at several triple points. The 
fundamental model with continuous spins and the Potts-type 
model with discrete spins share many properties, although in 
the absence of a centre symmetry they need not be in the same 
universality class. The continuous model exhibits a transi- 
tion from the symmetric to the ferromagnetic phase with the 
same behaviour of the Polyakov loop as in G2 gluodynam- 
ics, namely a steep jump from a small (but non-vanishing) 
Polyakov loop to a loop near the identity of G2. 

The classical, mean field and Monte-Carlo analysis all lead 
to a coherent and consistent picture for both 3-dimensional 
effective theories. In particular, the prediction of the mean 
field approximation for (^7) and (|x7,e — X7.o|) 1S m excel- 
lent agreement with the corresponding results obtained by de- 
tailed Monte-Carlo simulations. This parallels our findings 
for SU(3) in [18] and probably is due to the existence of tri- 
critical points which lower the upper critical dimension in the 
vicinity of these points. 

As concerning the relationship between the continuous ef- 
fective models to the underlying G2 gluodynamics we plan 
to apply inverse Monte-Carlo techniques, preferably with de- 
mon methods, to determine the dependence of the coupling 
in the fundamental model on the gauge coupling. We hope 
to present the resulting curve Xt((3), Ai4(/3) in a forthcoming 
publication. However we anticipate that for the critical Yang- 
Mills coupling /3 C this curve will cross the transition line be- 
tween the symmetric and ferromagnetic phases at small A14 
and negative A7 « —0.1 in Fig. 5. 
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